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Abstract. The Abelian sandpile growth model is a diffusion process for con- 
r" "H figurations of chips placed on vertices of the integer lattice % d , in which sites 

C ) with at least 2d chips topple, distributing 1 chip to each of their neighbors in 

O^l the lattice, until no more topplings are possible. From an initial configuration 

t I consisting of n chips placed at a single vertex, the rescaled stable configuration 

seems to converge to a particular fractal pattern as n — > oo. However, little has 

been proved about the appearance of the stable configurations. We use PDE 
^— ( techniques to prove that the rescaled stable configurations do indeed converge 

j-^^ to a unique limit as n — > oo. We characterize the limit as the Laplacian of the 

solution to an elliptic obstacle problem. 

j2 1- Introduction 

The Abelian sandpile growth model is a diffusion process for configurations of 
chips placed on vertices of the d-dimensional integer lattice 1 d . When a vertex 
in Z d topples, it loses 2d chips itself, adding one chip to each of its 2d neighbors 
in the lattice. For any finite nonnegative initial distribution of chips, one can 
topple vertices with at least 2d chips until no more such topplings are possible, 
i-H and, remarkably, neither the final configuration nor the number of topplings which 

occur at each vertex depend on the order in which topplings are performed — the 
process is in this sense Abelian, as first noticed in [5J[5]. 

Introduced in \U, the sandpile growth model has been the subject of extensive 
study over the past two decades, both on 7L d and on general graphs. Nevertheless, 
many of the earliest and most fundamental questions regarding the sandpile remain 
I unanswered (see e.g. |12| ). Some of the most nagging concern the final configuration 

r* s n : % d — > {0, 1, . . . , 2d— 1} obtained from starting with n chips placed at the origin 

. £^ of the integer lattice. Terminal configurations for several values of n and d = 2 are 

shown in Figure [l] Rescaled by a factor of n x l d , the piles seem to converge to a 
fractal-like limiting image, but essentially no progress has been made in proving 
precise things about this limit. (For example, there is no proof that there are 
regions in the limit which have constant value 3, no proof that its boundary does 
not converge to a circle, etc.) In this paper, we show at least that the limit exists 
as an object to study — the rescaled sandpile does converge. 

Theorem 1.1. The rescaled sandpiles s n (x) := s„(n 1 / d a;) converge weakly-* to a 
function s £ L°°(M. d ) as n — > oo. Moreover, the limit s satisfies Jl rf s dx = 1, 
< s < 2d - 1, and s = inR d \B R for some R>0. 
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FIGURE 1. Sandpiles started from n — 10 3 , n = 10 4 , n — 10 5 , 
and n = 10 6 chips placed at the origin of Z 2 , rescaled by factors 
of n 1 / 2 . Sites with 0, 1, 2, and 3 chips are represented by white, 
light gray, dark gray, and black, respectively. 



See Remark 4.3 for the characterization the limiting s e L°°(R ) as the Laplacian 
of the solution of an elliptic obstacle problem. 

We remark that weak-* convergence seems to be the correct notion of conver- 
gence for the rescaled sandpiles. Recall that L°°(M d ) denotes the space of bounded 
measurable functions on K d and that Co(K d ) denotes the continuous functions on 
R rf with compact support. A sequence of functions s n G L°°(Mr) converges weakly-* 
to a function s G L°°(R d ) in L°° (R d ) if 

s n <pdx-^ / sipdx as n —> oo, 

for all test functions <p G Co(IR !i ). If we look closely at the gray regions of the 
sandpiles (see Figure [2|, we see that they consist of rapidly oscillating patterns, 



as first noticed in 13 . It seems clear that no kind of pointwise convergence of 
the sandpiles can hold, assuming this behavior persists as n — > oo. Instead, the 
sandpiles converge in the sense that the oscillating regions converge to their average 
value; there is some limiting image s in which the "colors" of points are real num- 
bers, approximated in the sequence s n by patterns of integers. The convergence 
of a sequence of bounded functions in local average value is precisely what weak-* 
convergence in L°°(M <i ) is designed to capture [fj]. 
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Figure 2. Periodic structure in the "gray" regions of the n = 
16 ■ 10 6 sandpile on Z 2 . Notice the patterned triangular shapes. 



Note that the properties of s stated in the second part of the theorem are to 
demonstrate that the weak-* convergence is not for trivial reasons such as a poor 
choice of the rescaling: the mass of the limit is 1 and is not being "rescaled away" ; 
similarly, the rescaled sandpiles are contained within a ball of a fixed size. The best 
known upper bound on the diameter of the this ball [Tl] is that for every e > 
there is a C = C(e, d) > such that 

(1-1) { s n > 0} C B^ d _ e ^ Bl ^-i/d n l/d +c 

for all n > 0. Here \B\\ denotes the volume of the unit ball. 

Although our goal will be the characterization of final chip configurations, we 
will really be studying the number of topples which occur at each vertex as an 
initial configuration stabilizes. Given an initial configuration of chips and some 
finite sequence {pi} C T. d of points at which we topple piles of chips, the sequence 
is legal if the sequence never topples any vertex with fewer than 2d chips, and 
stabilizing if the final configuration has fewer than 2d chips at each vertex. Given 
any legal sequence pi , . . . , p s and any stabilizing sequence q% , . . . , qt both beginning 
from some fixed intial configuation, we must have p\ = for some i%, since there 
must be at least 2d chips at the point p\ in the initial configuration, since {p{\ is 
legal, and these chips must topple in any stabilizing sequence. Since the permuted 
sequence q^, q±, . . . , q^—x, q^+i, ■ ■ • , qt is also stabilizing, we can apply the same 
argument again from the configuration obtained after toppling just at p\ — q^ , to 
get that there must be a qi 2 — p2- Continuing in this manner, we obtain that the 
sequence {pi} is a permutation of some subsequence of {qi}. 

The above argument, a simplification of what appeared in jHJ, implies that for 
any fixed initial configuration of chips, any two legal stabilizing sequences must be 
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permutations of each other. Coupled with the fact that any finite initial configu- 
ration of chips on Z d will eventually stabilize after some legal sequence of topples, 
this implies that there is a well-defined odometer function v : Z d — > N associated to 
any initial configuration, which counts the number of topples which will occur at 
points of the lattice in any legal stabilizing sequence of topples. This also implies 
the Abelian property of the sandpile, as the final chip configuration s is uniquely 
determined by the odometer function. Note that it can be computed simply as 

(1.2) a(x)=ri(x) + ^2(v{y)-v{x)), 

where rj is the initial configuration of chips, and the sum is taken over the 2d 
neighbors y of x in the lattice Z d . Thus we have 

s = 7] + A 1 ?;, 

where A 1 is the (2d + l)-point discrete Laplacian on % d . 

The starting point of our proof is the least action principle formulated in |7j. 
Suppose that u is the odometer function for some initial configuration r] of chips, 
and v : Z d -> N is any function that satisfies 

X] + A 1 v <2d-l. 

The least action principle states that we must have u < v. Note that this follows 
immediately from the fact that legal sequences are permutations of subsequences of 
stabilizing sequences, since u corresponds to a stabilizing and legal sequence, while 
v corresponds to a stablizing, not necessarily legal sequence. 

The important consequence of the least action principle for us is that the odome- 
ter function associated with a configuration is the pointwise minimum of all func- 
tions v : Z d —> N satisfying n + A 1 v < 2d — 1, where n : Z d — s- N is the initial 
configuration of chips. In our case, if v n is the odometer function resulting from an 
initial configuration of n chips at the origin, we have 

(1.3) v n = mm{v : Z d -> N | n5 + A 1 v < 2d - 1}, 

where Sq is the characteristic function of the set {0} C Z d . The final configuration 
of chips is then 

s„ := n5 a + A 1 v n . 

This description of v n and s n , together with standard estimates for the (2d+ 1)- 
point Laplacian A 1 , will easily give convergence of the rescaled s n and the rescaled 
odometer function v n (x) := nT 2 l d v n (n}-l d x) along subsequences — > oo. 

To obtain convergence, we show that the limiting s and v are independent of the 
choice of subsequence. Assuming there are two distinct limits v and v', we use the 
regularity theory of the Laplacian to select a point x € M. d \ {0} where v and v' are 
sufficiently smooth and v — v' is strictly concave. We then select sufficiently close 
finite approximations v ~ v n and v' ~ v n i and use u„< to "lower" v n , contradicting 
the least action principle for v n . The difficultly lies in the fact that v n and v n > 
are defined on possibly incompatible lattices. We overcome this obstacle using an 



approximation argument (see Lemma 4.1 1 that allows us to change the scale of v n 
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2. Preliminaries 

2.1. Notation. We write h := n^ 1 /" 1 . We define the open ball B r {x) := {y E R d \ 
|y — ar| < r} for r > and x E R d and define B r := B r (0). We write dQ for the 
boundary of an open set O C R d . A function u E C(M. d ) is twice differentiable 
at a point x € M d if there is a vector Du(x) E K d and a symmetric d x d matrix 
D 2 u(x) 6 such that 

u(y) = + Dit(x) • (y - a;) + ^(y - xfD 2 u(x)(y - x) + o(\y - x\ 2 ). 

We write y ~/, x when y — x £ KL d := {ftx | a; € Z d } and |y — x| = ft,. Note that 
y ~h a? if and only if y = x ± Zie^ for some coordinate vector e,. Given a subset of 
the lattice £7 C /iZ d , we write 

d h E := {y E hZ d \E \y ~ hX for some x E E} 

for its lattice boundary. 

2.2. Interpolation convention. Throughout this article, we consider sequences 
of functions u n : hZ d — > M. We implicitly extend all such functions to M. d via 
nearest-neighbor interpolation. That is, we set 

u n (x) := u n (hlh~ 1 x~\) 

to be the value of u n at the lattice point h[h~ l x] E KL d closest to x E M. d . Here [•] 
denotes coordinate- wise rounding to the nearest integer (rounding down for ties). 
This convention allows us to make sense of statements like \ 4u6 C(IR d ) locally 
uniformly as n — >■ oo" when the functions u n are only defined on the lattice KL d . 

2.3. The Laplacian. The Laplacian 

d d 2 

Au := tra,cc(D"u) = > —^ru, 

and the (2d + l)-point discrete Laplacian 

A h u(x) :=i X>(y) -«(*)) 
y~h% 

play a central role in our analysis of the sandpile. Recall that if tp g C oa (M. d ), then 
A' 1 ^ — s- A(p locally uniformly in R d as h — >• 0. 



From 



10 , we know that the discrete Laplacian A h has a "fundamental solution'' 



$„ : KL d -> R that satisfies 

— n if x = 0, 



if x ^ 0, 



for all x E h!L d . The function $„(x) can be realized as a normalization of the 
expected number of visits of a random walk to the point x on the lattice hZ d , 
adjusting for recurrences in the case d = 2. As n — > oo, — > $ e C 00 (M d \ {0}) 
locally uniformly in M d \ {0}, where 



$(x) := 



(d(d - 2)\B 1 \)~ 1 \x\ 2 - d ifd>3, 
-(2 7 r)- 1 log|x| if d = 2, 
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is the "fundamental solution" of the Laplacian. We use the convergence of <!>„ — > $ 
to resolve the formation of a singularity in w„ as n — > oo. 

The discrete Laplacian is monotone: A h u(x) is decreasing in u{x) and increasing 
in u(y) for any y x. An obvious consequence of monotonicity is the following 
estimate for the discrete Laplacian of the point-wise minimum of two functions: 

Proposition 2.1. If u,v : hit — s- R, w := m.in{u,v}, and w(x) = u(x), then 
A h w(x) < A h u{x). 

A less obvious consequence of monotonicity is the fact that A h has a maximum 
principle [9j: 

Proposition 2.2. For u,v : hZ d 4i£C hZ d finite, and A h u > A h v in E, we 

have rn.Bx.Qh E (u — v) > maxs(ii — v). That is, the maximum difference must occur 
on the lattice boundary. 

Finally, we recall a consequence of the standard a priori estimates for the finite 
difference Laplacian and the Arzela-Ascoli theorem [9]: 

Proposition 2.3. If u n : KL — > K is a sequence of functions that satisfies 

max|u„| < C and max lA^Wnl < C, 

then for every sequence — > oo there is subsequence nk 6 — > oo and a function 
u G C(Bi) such that u Hk — > u locally uniformly in B\ as j — > oo. 

After we make some additional preparations, the above two propositions will 
allow us to show that the v n converge along subsequences. It remains to recall 
some standard facts about the continuum Laplacian which we use to study the 
regularity of the possible limits of the v n . 

If SI C M. d is open and s g L°°(S1), then we say that u € C(S1) is a weak solution 
(see [8]) of the continuum Laplace equation 

Am = s in SI, 

if and only if 

(2.1) / uAip dx = / sip dx 

Jn Jn 

for all test functions <p £ C^°(S1). Note that if u and s happen to be smooth, then 
we can we can intergrate by parts to show that this is equivalent to Au = s holding 
pointwise. 

The following is a consequence of a much stronger result from the theory of 
singular operators [3j Chapter III Theorem 4]: 

Proposition 2.4. If s £ L°°(Bi) and u £ C(Bi) is a weak solution of Au = s in 
Bi , then u is twice differ entiable and Au = s almost everywhere in Bi . 

Having bounded Laplacian also implies that a function must be strictly concave 
on a set of positive measure in any neighborhood of a strict local maximum [2j 
Theorem 3.2]: 

Proposition 2.5. If s £ L°°(Bi) and u £ C(Bi) is a weak solution of Au = s in 
B\, and sup Bl u > sup aBl u, then the set 

{x £ Bi | u(x) > supw, u is twice differentiate at x, and D 2 u(x) < 0}, 
has positive measure. 
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3. Convergence along subsequences 
Note from Section [I] and in particular ( |1.1| , that the sandpile s n satisfies 

(3.1) / s n dx — n, < s n < 2d — 1 and {s n > 0} C B h -i R , 

for some R > independent of n. These facts, together with the least action 
principle and the results in Section [2j are enough to prove convergence of the 
rescaled sandpile along subsequences. 

Recall that s n {x) := s(h~~ 1 x) and v n (x) := h 2 v n (h~ 1 x), and define 

Since A' 1 is linear and A <&„ = —ndo, we have 

A h w n = n<5 + A h v n = s n in M d . 

Lemma 3.1. For every sequence — > oo there is a subsequence rife. — > oo and 
functions w £ C{R d ) and s £ L°°(R d ) such that w is a weak solution of Aw = s in 
M. d , W nk . converges locally uniformly to w in C(M. d ), and s„ fc . converges weakly-* 
to s in L°°(R d ). 

Note that this lemma holds for any sequence w n : KL d — > M. such that w n and 
A h w n are locally uniformly bounded. Uniqueness of the limit is where our analysis 
is particular to the sandpile. 

Proof. Since {v n > 0} C Br and v n > 0, we therefore have 

w n > — <i>„ in KL d and w n = — $„ on hZ d \ B R . 
If we set E := hZ d n Br, then the test functions 

tp(x) := \x\ 2 — (R + h) 2 + inf -$„ and ip(x) := sup -$„, 

d h E d h E 

satisfy tp < w n < ip on d h E and 2c? = A h (p > A h w n > A h ip = in E. It follows 
from the maximum principle (Proposition 2.2 1 that 

\x\ 2 - (R + h) 2 + inf -$„ < w„(x) < sup -$ n for all x e B R . 

9 h E Qh E 

For any R! > R, we know that <!>„ — > $ uniformly in \ .Br as n — ^ oo. Thus 
the sequence w n satisfies 

\\Wu\\l^(B r ,) < C 

for some C = C(R') > and all n > 0. Since we also have lA^Wnl < 2d — 1, 
Proposition 2.3 implies the existence of a subsequence rife, and a function w £ C(R d ) 
such that w nj — > w locally uniformly as j — > oo. Since any uniformly bounded 
sequence of functions converges weakly-* along subsequences [6j Page 7], we may 
select a further subsequence and a function s £ L°°(M. d ) such that s nj converges 
weakly-* to s as j — > oo. 

For any tp £ Cfi°(R d ), we have 

(3.2) / tps nj dx= / tpA hi w nj dx = / (A h 'tp)w n . dx, 

where the second equality comes from the discrete integration by parts formula, 
which an be carried out in this case by writing A h ^ w n as a finite sum that commutes 
with the integral. 
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Since w n . — > w locally uniformly in R d , it converges uniformly on any closed 
ball, and therefore uniformly on any bounded neighborhood of the support of if. 
Since we also have A h ip — > Aip uniformly in we see that the right-hand side of 



(3.2) converges to J Rd wAipdx as j —> oo. The left-hand side converges to J Rd sipdx 
by the definition of weak-* convergence. Thus w is a weak solution of Aw = s in 
R d by definition. □ 



An immediate consequence of Lemma 3.1 is that v n . — ¥ v := w + $ locally 
uniformly in IR'' \ {0} as j — > oo. However, what we obtained is strictly stronger, 
since it allows us to resolve the structure of the singularity of v. 

4. Convergence 

As discussed in Section [TJ we would like to compare two limits v nk — > v and 
v n i — » v' along different sequences rih — > oo and n' k — > oo by comparing their finite 
approximations, in spite of the fact that their approximations v nk and V n ' may 
not be defined on the same scale. To deal with this, we could try construct an 
asymptotic expansion of the convergence v rik — ¥ v in a neighborhood of each point 
x £M. d \ {0} where v is twice differentiable. That is, we could try to find a radius 
r > and a function u : 7L d — > Z such that A x u < 2d — 1 and 

(4.1) v nk (y) = v(x) + Dv(x) ■ (y - x) + hlu(h^(y - x)) + o(\y - x\ 2 ) + o(l), 

for all y E B r (x) and k > 0. (Note that o(l) is with respect to k — > oo while 
o(\y — x\ 2 ) is with respect to the distance \y — x\ — > 0.) The existence of such an 
expansion would make it easy to approximate v on different scales in B r (x). Indeed, 
while v only has approximations v 1lk defined on the scales hk, the right-hand side 



of (4.1 1 is valid for any scale h > 0. 



Unfortunately, we do not know how to prove that an asymptotic expansion of 



the form (4.1) exists in general. Instead we construct a one-sided expansion, which 
is sufficient for our purposes. Note that this lemma is applied below to v n that 
have been translated and shifted and in particular may take negative values. 

Lemma 4.1. Suppose v n : Z rf — > Z is any sequence of functions such that A 1 v n < 
K G Z. Suppose v G C(B r (xo)) for some r > and the rescalings v n (x) :— 
h 2 v n (h~ 1 x) converge uniformly v nk — > v in B r (xo) along some sequence — > oo. 
If v is twice differentiable at Xq, then for every e > there is a function u : Z d — > Z 
such that 

A 1 u(x)<K and u(x) > ^x t (D 2 v(x ) - el)x forallxeZ d . 

Proof. Replacing v n if necessary by 

v' n (x) := v n (x + L^ _1 a:ol ) - v n (lh~ 1 x Q ~} ) - [h^ 1 Dv(x )~\ ■ x, 

we may assume that xo — 0, v(0) — 0, Dv(0) = 0, and v n (0) = 0. 

Since v is twice differentiable at 0, we can make r > smaller and select a large 
n = rife such that 

sup \v n - ip\ < er 2 , 

B r 

where <p(x) := ^x t D 2 v(0)x. Noting that h 2 (p(h~ 1 x) — ip(x), we can undo the 
scaling of v n , obtaining 

sup \v n — ip\ < sh~ 2 r 2 . 
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Thus, if we define tp(x) := ^x (D v(0) — 32e/)x, we have ip — ip = 16e|a;| and 
therefore 



(4.2) 



in Bh-i r /4, 



\v n < ip + 2eh~ 2 r 2 

\v„ > ip + 3eh~ 2 r 2 in B h -i r \ B h -i r/2 . 

We define u as an overlapping pointwise minimum of translated and tilted copies 
of v n , 

u(x) := mm{v rit y(x) \ y £ 7L d n B h -i r (x)}, 

where 

v n ,y{x) := v n (x - y) + [Dip(y)] ■ (x - y) + [ip(y)~\. 

The inequalities (4.2) guarantee that the overlapping works out correctly. Indeed, 
for x £ Bf l --L r i i {y), we compute 

v n ,y(x) =v n (x-y)+ [Dip(y)] ■ (x - y) + [ip(y)] 

< i/j{x -y)+ L^(j/)l -(x-y)+ LV>(y)l + 2eh~ 2 r 2 

< ip(x ~y) + Dip{y) ■ (x - y) + tp(y) + 2Eh~ 2 r 2 + d lj2 hr x r + 1 
= ip(x) + 2ehr 2 r 2 + d 1/2 ^V + f . 

A similar computation shows 

v n , y > + iebT 2 r 2 - d x l 2 hT x r - f in B h -i r (y) \ B h -i r/2 (y). 
By making n larger, we may assume that h~ 1 r/4 > 1 and 



2eh- 2 r 2 + d 1/2 h- 1 r + 1 < 3e/i~ V 



d x ' 2 h- 



1. 



It follows that 
and therefore 



u(x) = min{i;„^(a;) | y £ 7L d n B h -i r/2 {x)} 1 
min{v n ^ y | y £ Z d f] B 3h -i r/4: (x)} in B h ~i r/4 (x), 



for all x £ 7L d . Since A 1 ^^ < K and h 1 r/4 > 1, Proposition 2.1 implies that 
A 1 u(x) < K. Moreover, it is easy to check that u — ip is bounded from below. □ 

We are now in a position to prove Theorem |1.1| In fact, we prove something 
slightly stronger. 

Theorem 4.2. There are functions w £ C(R d ) and s £ L°°(R d ) such that s n 
converges weakly-* to s in L°°(R d ), w n converges locally uniformly to w in C(R d ), 
and w is a weak solution of Aw — s in Mr. Moreover, the function s satisfies 
J md sdx = 1, < s < 2d — 1, and s = in IR d \ Br for some R > 0. 



Proof. Since Lemma 3.1 gives existence of limits along subsequences, we need only 
prove uniqueness of the limiting s and w. By Proposition 2.4 Aw = s almost 
everywhere, so it is enough to show that w is unique. Suppose w,w' £ C(M. d ) are 
distinct and that w nk — > w and w n ' k — > w' along subsequences locally uniformly in 
E d as k — > oo. Since w = w' = — $ outside of Br for some R > 0, we may assume 
without loss of generality that 

sup(w — w') > = sup(ui — w'). 
b r ob r 
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According to Proposition 2.5 we may select a point a E R d such that a^O, w(a) > 
w'(a), both w and w' are twice differentiable at a, and D 2 w(a) < D 2 w'(a) — 2el 
for some e > 0. 

Define v := w + <f> and v' := to' + Since n / 0, ?;„' — > i/ uniformly in a 



neighborhood of a as /c — >■ oo and v' is twice differentiable at a. Using Lemma 4.1 
we may select a u : Z d -> Z such that 

A 1 u(x) <2d-l and u(x) > -x t {D 2 v'{a) - f?i> for all a: G Z d . 

Since D 2 v'(a) - el > D 2 v(a) + el it follows that 

> -x t (D 2 v(a) + el)x for all x E Z d . 

We use a shifted version of u to "lower" u njfc . For the a chosen above, we define 

h k := ra^ 1 , a fc := [/i* 1 ^, 

and 



(a; - a k ) + [h, Dv(a)] ■ (x - a k ) 



ti\~ . ,, , ......... : . , 5 

«,,,(V, :=<| -u(0)+v nk (a k )-l ^'VW' 

1 + max Z d w„ fe otherwise. 

We claim that for small r > and large n k , the function u := mm{v nk ,u nk } 
contradicts the least action principle for v nk . Since 

v{a k ) = u nk (a k ) = v 7lk {a k ) - 1, 

it is enough to show that uSq + A 1 ^ < 2d — 1 and v > 0. 

To show that n<5o + A 1 ?; < 2d — 1, it is enough, by Proposition 2.1 to show that 
v Uk < u n . in B.-i (ofc) \ B,-i /2 (a k ) and ^ B,-i ,Ja k ) for small r > and 

fc k I k I 

large rife. Since w is twice differentiable at a, we have, for small r > and large n^, 
w nfc (a:) < v nk (a) + Dv(a) ■ (x - a) + -(x - a) t D 2 v(a)(x - a) + ^r 2 , 

for all x E B r (a). Undoing the scaling, we discover that 

1 e 

Vn k {x) < v nk (a k ) + h~ 1 Dv(a) ■ (x - a k ) + -{x - a k ) t D 2 v{a){x - a k ) + -^h k 2 r 2 , 

for all x E B,-i (a k ). Since 

k 

1 e 
u(x) > -x t D 2 v(a)x + -^K 2 r 2 , 

for all x E B h -i r \ B h -i r ^ 2 , we see that 

u nk (x) > v nk (x) + £ -h- k 2 r 2 - d^h^r - u(0) - 1, 
for all y E B,~i (a k ) \ B,-i ,Ja k ). Since the error term on the right-hand side is 

k k I 

positive for large n k , we see that 

u nk > v nk in B h -i r (a k ) \ B, -i , 2 (a k ), 

k k l 

for small r > and large n^. 
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To show that v > 0, it is enough to show that u nk > in B h -i (ak) for small 

k 

r > and large Ufc. Note that since v is twice differcntiable at a, — > t> in B r (a), 
and w(a) > u>'(a), we have 



w„ t (a) +Dv(a) ■ {x - a) + -{x - af D 2 v{a)(x ~ a) > > — ^ > 0, 
for small r > 0, large n^, and x G B r (a). Thus, for x G B h -i r (a k ), we may compute 
u„ fc (x) = u(x - a*,) + [fefe 1 -D«(ffl)l • (a: - flfc) - u(0) + w„ fc (a fe ) - 1 

> i(a; - a k ) t D 2 v{a){x ~ a k ) + h^ 1 Dv(a) ■ (x - a k ) + v nk {a k ) 

-u(0) - 1 -d 1/2 \x~a k \ 

>K^-u{Q)-l-h- k ^ 

> 0, 

for small r > and large n k . 

The remaining assertions about s follow immediately from (3.1 1 and our choice 
of scaling. □ 

Remark 4.3. For those readers familiar with viscosity solutions of fully nonlinear 
elliptic equations (see (2|), we point out that the above proof shows that the limit 
w n — > w is the unique solution of the obstacle problem 

(4.3) w := infjV <E C*(K d ) | w' > -$ and Aw' < 2d - 1 in R d , 

and G{D 2 w' + D 2 <S>) < in R d \ {0}}, 

where 

G(A) := inf{s G K | there is a u : Z d -> Z such that for all y E Z d 

A 1 u(y) < 2d - 1 and 11(2/) > iy*(,4 - s/)y}, 

and we interpret the differential inequalities in the sense of viscosity. In fact, our 
proof that the limit w is unique is exactly a proof that w is the unique solution of 
(4.3) with the viscosity solution terminology stripped out. 

Without too much difficulty, one can show that the operator G is Lipschitz and 
elliptic. However, we have no reason to believe that G is uniformly elliptic, and 
thus the existing theory of fully nonlinear elliptic obstacle problems does not apply 
directly. Whether this theory can be re-worked for G and whether the characteri- 
zation (4.3) is useful remains to be seen. 
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